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The problem of the figure of the earth has, as far as I know, only 
received one solution, namely, that of Laplace.* His solution involves 
an hypothesis as to the law of compressibility of the matter forming 
the planet, and a solution involving another law of compressibility 
seems of some interest, even although the results are not perhaps so 
conformable to the observed facts with regard to the earth as those of 
Laplace, f 

The solution offered below was arrived at by an inverse method, 
namely, by the assumption of a form for the law of the internal 
density of the planet, and the subsequent determination of the law of 
compressibility. One case of the solution gives us constant compres¬ 
sibility, and another gives the case where the modulus of compressibility 
varies as the density, as with gas. 

It would be easy to fabricate any number of distributions of 
density, any one of which would lead to a law of compressibility 
equally probable with that of Laplace; but the solution of Clairaut’s- 
equation for the ellipticity of the internal strata of equal density 
seems in most cases very difficult. Indeed, it is probable that Laplace 
formulated his law because it made the equation in question integrable,. 
and because it was not improbable from a physical point of view. 

The following notation will be adopted :— 

For an internal stratum of equal density let— 

r be the radius vector of any point, 
a the mean radius of the stratum, 
e the ellipticity, 
w the density, 

0 colatitude from the axis of rotation, 
jp the hydrostatic pressure at the point r, 0. 

For the surface let tr, a, t, It) denote the similar things. 

Let M be the mass of the planet, p its mean density, iv the angular 
velocity of rotation, m the ratio w z /^7rp. 

# Since this paper was presented I haye seen a reference to a paper by the late 
M. Edouard Roche, in vol. i of the Memoirs of the Academy of Montpellier (1848), 
in which the problem is solved, when the rate of increase of the density varies as 
the square of the radius. See Tisserand, “ Comptes Rendus,” 23rd April, 1883. 

f Laplace’s hypothetical law of compressibility arises from a law of internal, 
density for which the problem had previously been worked out, as an example, by 
Legendre. See Todhunter’s “History of the Figure of the Earth,” vol. ii, pp. 117 
and 337. 
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Let h be tbe ratio of tbe density of tbe stratum a to tbe mean 
density of all tbe matter sitnated inside tbat stratum, and fe tbe 
surface value of lcJ /? 

Let C, A be tbe greatest and least principal moments of inertia of 
tbe planet about axes tbrougb its centre of inertia. 

Let ce be tbe ellipticity which tbe surface would bave if tbe planet 
were homogeneous witb density p, so tbat 

Tbe condition tbat tbe surface of tbe planet is a level surface is 
satisfied by— 

C-A=-|Ma 3 (e-im) ...... (1). 

Tbe condition tbat tbe internal surfaces are also surfaces of equi¬ 
librium demands tbat e should satisfy Clairaut’s equation—- 

(S- 6 y)f>“ !,i “ +2 “ ! (l + ;) =0 ■ ■ • • $'>■ 

\tvtv Lv / J U \U/W/ tv/ 


It may be proved from (2), and tbe consideration tbat iv must 
diminish as a increases, tbat e cannot bave a maximum or minimum 
value. 

Also it may be shown tbat tbe constants introduced in the integral 
of this equation must be such tbat— 


t 4 e 2sa da 


(3). 


when a is put equal to a after differentiation. 

Tbe mean density is given by— 

a 3 p=3j* wo?da .. (4). 

And 6=-. 

P 

Neglecting tbe ellipticity of tbe strata, we bave tbe moment of 


inertia about any diameter of tbe planet given by— 

C ~ 3 7r j' ‘wvft&u .. (5). 

Tbe ratio of (1) to (5) gives tbe precessional constant. Tbe pres¬ 
sure and density are connected by tbe equation— 

j* (-^E-\-47rwa\da-\-^—\ wa?da~ 0 .... (6). 

J a\w da J a Jo 


* 7c, It, are the reciprocals of f, f, according to the notation adopted in Thomson 
and Tait’s “b T at. Phil.” (edit, of 1883), § 824. 
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Now if m be a function such that wdjp—dm, the differentiation of 
(6) leads to— 

d ^ + ^.\\va>da =0.(7). 

da a z J o 

and a second differentiation to— 

(ma) +47twu=0 ....... (8). 

da 2 


It is well known that Laplace assumed that the modulus of com¬ 
pressibility of rock varies as the square of the density. Since this 
modulus is wdjp/dw, Laplace’s hypothesis makes m proportional to w, 
and the equation (8) is at once soluble. 

After the determination, of w as a function of a, the solution of all 
the other equations follows. 

In this paper I propose to find a new solution, and to compare the 
results with those of Laplace. 

In order to simplify the analysis let the unit of length be equal 
to the mean radius a of the planet, and the unit of time be such that 
the surface density iu of the plant is also unity. 

Now let us assume that the law of internal density is— 


w=a~ n .. (9). 

Then the mean density of all the matter lying inside of the stratum 
a is a~~ n l (1—^n). Hence, by definition we have— 

........ ( 10 ). 


Thus we see that h is a constant for all strata, and therefore also 
for the surface. In Laplace’s theory k is variable. With our assumed 
law of density and the special units, /> the mean density is equal to 
the reciprocal of k. 

It is clear that n must be positive, otherwise heavier strata lie 
above lighter, and it must be less than 3 in order to avoid infinite 
mass at the centre of the planet. 

Now let us find the law connecting pressure and density, and the 
modulus of compressibility. 

Equation (7) becomes 

dm . a l ~ u 

- = -f 7T--—, 

da Y—^n 

and by definition of ci and the assumption (9), 

dp __ 4,^ aY~ %l 

da 3 1 
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Integrating this, with the condition that the pressure vanishes at 
the surface, we have, 

P~‘TT\ rv7i ® 2<1— “’J) 


3(1—w)(l- 

27r 




3(1-«)(!-» 


[1— w~ 2 » ], 


whence the modulus of compressibility is 
dp d 1 

or: J —- ** cr _ 4 /} 


dw 


n (\—k n ) 

The case of n= 1 is interesting; it gives a constant modulus of 
compressibility equal to 27r, and the law of pressure p = 27r log w. 

If w be less than unity the compressibility, or reciprocal of the 
modulus, increases with the density, which is of course physically 
improbable. If n be greater than unity and less than 3, the com¬ 
pressibility becomes less the greater the density. The assumed law 
probably does not give such good results as those of Laplace, because 
the decrease of compressibility with increasing density is not 
sufficiently rapid. 

The range of n =3 to n —1 gives the results which possess most 
physical interest. 

In comparing results with those of Laplace there will be occasion 
to express the modulus as a length; that is to say, we are to find the 
length of a column of unit section whose weight (referred to the 
surface gravity of the planet) is equal to the force specified in the 
modulus. 

Now if g be gravity 

g=±7rp— -fT T- - . 


a ^ ^ 

Hence the modulus is =g afo X 


, 1 /w\ 3 ir 

n Vw ? 


the units a, tu being 


reintroduced to give the expression the proper dimensions. Now 
is a pressure, and therefore the length of the modulus is 

/ A 2n—1 

- y? j n . Thus the surface matter has a length modulus equal to a jn. 

Now let us find the ellipticity of the internal strata. 

Substituting for w from (9) in (2), we have, 

af± + 2(3-n)^--2,n-=0. 

dor da a 

If the solution be assumed of the form e=ca£, /3 must satisfy 
/3(/3—1) +2(3— ri)j3 —2^=0, 
/ 3=-(|- W )± v /{(4)2- w (3- W )}. 


whence 
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Now n(3—n) is a maximum when it is equal to (f) 2 , and therefore 
the square root can never become imaginary. 

From the sign of the last term in the equation for (3, it is clear that 
one of the values of ft is negative. Hence to avoid infinite ellipticity 
at the centre, the c corresponding to the negative root must be zero. 
Hence the solution of Clairaut’s equation (2) is 

0 ~ %oJ l ~ 5+ ^ [(f) 3 — n(3—n)]" 

The surface value of is clearly n— \/{(f) 3 —-*w(3—w)}.. 

z da 

Thus from (3) we have 


Then substituting for w from (9) in (5), 

o=- 8,r 


3(5—w) 


And since M=f7r^=f7r-—we have for the precessional constant 
1 s n 

0-A = |^ (e _ W . 

Now let us collect these results, and express them in terms of k 
instead of n. The solution is 

And the mass inside of any radius a is 

ole 

4-6 k 

W 3-3 k —1 

p=b 


div 


Aj(2~3 Jc) 

4-6 k 
q/j'i —3 k 


Sk(l-k) 

And when &=§,_£> = 27r log iv, w$£-~2i r, w = a~ J . 

dw 

The length of the modulus at the surface is 1/3(1— k) of the 
planet’s radius. 

e=zta h - sk+s 


——i—-ik + f v /{(f) 2 —k( 1 --k}„ 

e 

C= . ~ . §Ma 2 . 


2+3& 
C-A_2+3fc 


C 


3& 


0 -\m). 
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Any value from unity to an infinitely small value may be assigned 
to h, that is to say, we may have any arrangement of density from 
homogeneity to infinitely small surface density, but if h be greater 
than -J the compressibility increases with the density, which is 
physically improbable. 

The infinite density and infinite pressure, which occur in this 
solution actually at the centre, may be avoided by imagining the 
centre occupied by a rigid ;;spherical homogeneous nucleus, of very 
small radius and of density l//^u 3(1 “^. 

We have to compare this solution with Laplace’s. 

For this case k is not constant, and its surface value is fc. 

Let $=-, where k is a constant, being the arbitrary constant intro- 

K. 

duced in this solution ; and let 9 be the surface value of 
The solution is— 

6 sin# 

iv—~ —- 

$ sm 9 

And the mass inside of any radius u is — a Sk . 


jp=27ric i (w® — 1 ). 
W < ^=47ne 3 W 2 . 

dw 

h=z _ 

3(1 — # cot #) 


The length of the modulus at the surface is 1/(1 — 9 cot 9) or 3 
of the planet’s radius. 




0*1—k 

Vi-f 


ce_ e * 

t _ 6(1—fe) 


#6. 


C=[l —6(1—fe)0~ 2 ]§Ma 3 . 


C-A 

0 


1 

1—6(1—fe)0 -2 


(e—| m). 


The following table gives the numerical values of the solution, 
together with columns for comparison with the results of Laplace’s 
theory, for various values of the ratio of surface to mean density. 


# See Thomson and Tait’s <( Nat, Phil.,” 1883, § 824. 
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The case of k=0 gives the planet infinite mass at the centre, and 
the values are only inserted in order to complete the series. 


k. 

Length mod. 
at surface 
in terms of 

K as unity. 

e^Za x 

where 

X — 

ce 

F 

Laplace 

re 

Z 

C-A 

Laplace 

C-A 

po£ {iVy^r 1). 

where 

C (Z — &n). 

C(£ — \m). 

1*0 

00 

0 

000 

1*000 


•00 

1*667 

1 ‘67 

— 00 

*9 

3*333 

0 

132 

1 *066 

l 

’04 

1 '741 

1-71 

-4*667 

*8 

1 *667 

0 

293 

1*147 

l 

* 09 

1 *833 

1 ’77 

-1 *333 

*7 

1*111 

0 

488 

1*244 

7 

* 15 

1 *952 

1 ’83 

-0*067 

*6 

1-4-1 *2 

0 

722 

1*361 

1 

‘21 

2*111 

1’90 

0*333 

*5 

14-1*5 

1 

000 

1*500 

1 


2*333 

1 ’98 

0*667 

*4 

1-4-1 *8 

1 

322 

1 *661 

1 

’34 

2*667 

2 ’07 

0*889 

*3 

1-4-2-1 

1 

688 

1 *844 

1 

’41 

3*222 

2'17 1 

1*048 

*2 

14-2*4 

2 

093 

2 *047 

1 

’48 

4*333 

1 0 ’28 ! 

1*167 

*1 

14-2*7 

2 

532 

2 *266 

1 

’56 

7-667 

2’40 ! 

1*259 

*0 

14-3*0 

3 

000 

2*500 

1 

’65 

00 

2 '55 

1*333 

•667 

1*000 

0 

562 

1*281 

1 

’17 

2 *000 

1 '85 

I poc log w 

*333 

1-4-2 *0 

1 

562 

1 *781 

1 

’39 

3*000 

2 '13 

j 1*000 


Note .—The values in the two columns applicable to Laplace’s theory were found 
by graphical interpolation from a series of values given in “ Month. Wot.,” W.A.S., 
Dec., 1876, or Thomson and Tait, “Wat. Phil.” (1883), § 824'. 

In Laplace’s theory poc (w 2 — 1), and the modulus of compressibility oc iv l . In 
the present theory the modulus oc wv. 

The value it = *667 corresponds to constant compressibility, and it = *333 to 
gaseous compressibility. 


One of the grounds on which Laplace’s solution is held to be satis¬ 
factory is that if we take the value of ce , as determined by the known 
angular velocity and mean density of the earth, and the value of z as 
determined by geodesy, and find the value of k, the ratio of surface to 
mean density, which corresponds with the ratio ce/z, this same value 
of k is found to give a proper value to the coefficient of z —so as 
to obtain the observed precessional constant. To be more precise, m 
is found to be 1/289*66, which gives ce — 1/231*7, and z has been found 
to be approximately 1/295. These give ce/z~ 1*273, and this corre¬ 
sponds with k=l/2*06 —*49. This value of k, with the same values 
of z and m, gives the precessional constant as *0033, and Leverrier 
and Serret give its value as *00327. 

Now it appears remarkable that almost as good a correspondence 
is obtainable from my solution. The value ee/e—1*273 corresponds 
with k=*675, and when k = *675 the coefficient of z — -|m in the pre¬ 
cessional constant is 1*99, which gives the same precessional constant 
*0033. 

This value of k corresponds very nearly with constant modulus of 
compressibility, and with pressure determined by p — 27r log w. 
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It is claimed in favour of the Laplacian hypothesis that it corre¬ 
sponds to a surface density which is nearly a half of the mean 
density of the earth, and that we know that average rock has a den¬ 
sity of about 2*8. Also it is pointed out in Thomson and Tait’s 
“ Natural Philosophy ” that the length modulus of compressibility 
of the surface rock is about 1/4*4 of the earth’s radius, which is 
very nearly the observed length modulus of iron. 

These conditions are not well satisfied by the present solution, for 
the surface density is found to be *675, or 1/1*48 of the mean density 
of the planet, whence the specific gravity at the surface is 3*7; and the 
length modulus at the surface is equal to the planet’s radius. It is fco 
be admitted that this density is large, and that the substance is also 
highly incompressible. 

Thus in these respects the Laplacian hypothesis has the advantage. 
It seems to me, however, that too much stress should not be laid on 
these arguments. We know nothing of the materials of the earth,, 
excepting for a mile or two in thickness from the surface, hence it is 
not safe to argue confidently as to the degree of compressibility of 
the interior. T here seems reason to believe that there is a deficiency 
in density under great mountain ranges, and this would agree with 
the hypothesis that our continents are a mere intumescence of the 
surface layers. 

According to this view we might expect to find a rather sudden 
change in density within a few miles of the surface. Now in any 
theory of the earth’s density such a sudden change in the thin shell 
on the surface could not be taken into account, and the numerical 
value for the surface density should be taken from below the intu- 
meseent layer if it exists. Hence it is not unreasonable to say that a 
solution of the problem, which gives a higher surface density than 
that of rock, lies near the truth. I do not maintain that my solution 
is as likely as that of Laplace, but it is not to be condemned at once 
because it does not satisfy these conditions as to the density and com¬ 
pressibility of rock. 

The two cases which are given at the foot of the above table each 
possess an interest, the first of constant compressibility, because it 
corresponds with the case of the earth, and the second of modulus of 
compressibility varying as the density, because this is the gaseous 
law. 

With constant compressibility the internal ellipticity varies as the 
562 power of a , or nearly as the square root of the radius; with 
gaseous compressibility it varies as the 1*562 power of a, or nearly as 
the square root of the cube as the radius. 

A n um erical comparison of the case of constant compressibility 
with Laplace’s solution for gives the following results :—• 
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a — 

0 

ia. 

i a. 

•| a. 

a. 

(Laplace) ~ = 

0 

*812 

•844 

•902 

1-000 

(Constant e 
compress.) ^ ~ 

0 

•459 

*693 

*851 

1*000 


Thus the Laplacian solution attributes much higher ellipticity to 
the internal strata. The solution with constant compressibility in 
fact gives so large a proportion of the mass in the central region, that 
attraction has a greater influence compared with rotation, than in the 
solution of Laplace. 

[P.S. —If, as is not improbable, the increase of density in the interior 
of the earth is due rather to the heavier materials falling down to the 
centre than to great pressure compressing the material until it has a 
high density, then the determination of a modulus of compressibility 
would be fallacious, and it would be more logical to leave the 
expressions for the pressure and the density both as functions of the 
radius, without proceeding to eliminate the radius and to form an 
expression for the modulus of compressibility. I owe this suggestion 
to a conversation with Sir William Thomson.—December 19, 1883.] 




